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O . In the context of neutrino scenarios characterized by four (three active plus 

one sterile) neutrino species and by mass spectra with two separated doublets, 
we analyze solutions to the atmospheric neutrino anomaly which smoothly in- 
terpolate between v^ -^ v-j- and f^ -^ Vs oscillations. We show that, although 
the Super-Kamiokande data disfavor the pure v^ -^ Vg channel, they do not 
exclude its occurrence, with sizable amplitude, in addition to the v^ — > Vr 
Q> \ channel. High energy muon data appear to be crucial in assessing the relative 

^ \ amplitude of active and sterile neutrino oscillations. It is also qualitatively 

Q"^ ' shown that such atmospheric v solutions are compatible with analogous so- 

^P ■ lutions to the solar neutrino problem, which involve oscillations of Ve in both 

sterile and active states. 
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I. INTRODUCTION 

Four-neutrino (Au) models, involving the three known active states {ue, z/^, i^r) plus one 
hypothetical sterile state {ug), are being intensively studied in the context of the current 
neutrino phenomenology, since they can accommodate (through three independent mass 
square differences) the three sources of evidence for z/ flavor oscillations coming from atmo- 
spheric, solar, and accelerator neutrino experiments. In particular, 4z/ spectra with mass 
eigenstates organized in two doublets (2-|-2 models) seem to be favored by world neutrino 
data [H]. Updated reviews on such Au models and on their phenomenology can be found, 
e.g., in |0-§], to which the reader is referred for an extended bibliography and for details 
that are not repeated in this work. 

According to the conventional wisdom, 2+2 models are often assumed to imply, for the 
solar and atmospheric oscillation channels, either 

(la) 
(lb) 

or 

(2a) 
(2b) 

in addition to inter-doublet z/^ -^ z/g oscillations with small amplitude, required to explain 
the Liquid Scintillator Neutrino Detector (LSND) signal P,|^. 

Such simplifying assumptions are challenged by the recent Super-Kamiokande (SK) ob- 
servations, which tend to disfavor (dominant) oscillations into Us of both atmospheric muon 
neutrinos (at ~ 99% C.L. 0,|]) and solar electron neutrinos (at ~ 95% C.L., in combination 
with world solar z/ data [Q). The MACRO experiment also disfavors atmospheric u^ —>■ Ug 
oscillations (at ~ 95% C.L. [|1^). However, it seems premature to rule out the sterile neu- 
trino hypothesis on the basis of the present information only [11,1^: on the one hand, past 



z//.- 


->■ Ur (atmospheric 


l^e - 


-> Us (solar) , 


^M- 


-^ Ug (atmospheric 


Z^e - 


->■ Vr (solar) , 



experience on global data fits has shown that, sometimes, hypotheses rejected (accepted) at 
95% or 99% C.L. can surprisingly revive (die); on the other hand, the underlying 4z/ oscilla- 
tion pattern might be more comphcated [p|, p!3| , p!4| than assumed in Eqs. ( p!a] , p!b| ) or (pa| , |2b|) . 



In particular, instead of having decoupled active and sterile oscillation channels, one might 
have mixed (active+sterile) flavor transitions of the kind 



^-'^l -^ ^+ (atmospheric) , (3a) 

Ue -^ z/_ (solar) , (3b) 

where the states z/±, as discussed later in more detail, represent linear (orthogonal) combi- 
nations of Uj. and z/^ through a mixing angle ^, 

Z/+ = + cos C.h'r + sin ^ Ug , (4a) 

z/_ = — sin ^ z/^ + cos ^ Ug . (4b) 



The oscillation modes (|3a|j3BD represent generalizations of both modes (|T^,|1^) and (pa|,|2' 



to which they reduce for sin^ = and 1, respectively. For generic values of sin^, the final 



states in atmospheric u^^ and solar z/g flavor transitions are linear combinations of Us and z/,-, 
with coefficients to be constrained by experiments. 

In this work, we quantitatively study atmospheric four-neutrino oscillations in the context 
of 2+2 spectra, for unconstrained values of sin^ = (z/_|_|i/s). It is shown that the state u^ 
(into which i/^ oscillates) can have a sizable Us component, which is testable through high- 
energy atmospheric muon data. Our atmospheric z/ results are also qualitatively compared 
with recent 4z/ solutions to the solar neutrino problem |TB|,|T^, which are compatible with a 
large Ug component of i'_. It is shown that both solar and atmospheric neutrino oscillations 
are consistent with sizable flavor transitions to a sterile state. 

The structure of this paper is as follows. In Sec. II we introduce the theoretical 4z/ frame- 
work, and in Sec. Ill we show the corresponding graphical representation of the parameter 
space. Section IV reports a preliminary discussion of the SK experimental information, 
while Sec. V is devoted to a thorough SK atmospheric ly data analysis, from which we derive 
constraints on Ai/ mass-mixing parameters. Their consistency with solar neutrino data is 
discussed in Sec. VI. Remarks on alternative 4i> scenarios are made in Sec. VII. Finally, 
conclusions and perspectives are reported in Sec. VIII. 



While this work was being completed, our attention was brought to the paper |[16 
where 4z/ solutions to the atmospheric neutrino anomaly have also been worked out. The 
4z/ scenario considered in |^ contains two additional parameters as compared with ours (a 



mixing angle S23 and a CP violation phase 61 [0]), and most of the results shown refer to 
such parameters, so that a direct comparison with our results is difficult. However, we agree 
with [|16| on some general features, namely, that atmospheric z/^ can have a large transition 



amplitude to Ug (in addition to z/^ -^ Ur), and that such transitions are consistent with solar 
neutrino data. 



II. THEORETICAL FRAMEWORK 

In this section we define the theoretical 4z/ framework and the notation used in our 
analysis. We also motivate a few data-driven assumptions, which considerably simplify both 
the analysis and the understanding of the solar and atmospheric z/ oscillations. 

A. General conventions 

We order the neutrino flavor and mass eigenstates in column vectors as 

Ua = {l^e, l^s, I'^i, I'tY' , (5) 

Vi = (z/i,Z/2,Z/3,Z/4)^ , (6) 

respectively (T = transpose). Such vectors are related by a unitary mixing matrix U, 

Va = UaiVi . (7) 

The matrix of squared neutrino masses (in the z/j basis) is defined as 

M = dia.g{m\,ml,ml,m\) + \^1 , (8) 



where the term A^l reminds that the overall mass scale A is unconstrained, as far as os- 
cillation experiments are concerned. We define three independent mass square differences, 
phenomenologically related to the evidence for oscillations coming from solar, atmospheric, 
and LSND experiments, 

6m'^ = m\ — m\ (solar) , (9a) 

rn^ = ml — ml (atmos.) , (9b) 

M^ = ml- ml (LSND) . (9c) 

With the above definitions, the neutrino evolution equation in the fiavor basis, 

d 
«— Z/q = 'Hva , (10) 

is governed by the Hamiltonian 

H=^f/A^f/1" + V2G'^diag(iV„iiV„,0,0) , (11) 

where the first (kinematical) term [|r^ depends on the neutrino energy E, and the second 
(dynamical) term [0 depends on the neutrino weak interactions with background electrons 
(with density iVg) and neutrons (with density Nn) in matter. 

B. Approximations concerning Av masses and mixing 

Motivated by the current Au oscillation phenomenology (see f^-Q and references therein), 
we assume that the mass eigenstates are organized in two doublets (2+2 spectrum), 

mi ~ m,2 < m,3 ~ m,^ . (12) 

with three widely different mass squared differences, 

5m^ <t: m^ <t: M'^ . (13) 

Figure 1 shows our reference mass spectrum, with a lower "solar neutrino doublet" 
(z/i,z/2), and an upper "atmospheric neutrino doublet" {1^3,1^4), separated by a relatively 
large LSND mass gap. Other phenomenologically allowed 2+2 spectra can be obtained 
by interchanging the two doublets (M^ -^ — M^), as well as the two states in a doublet 
((5r7i^ -^ —5im?^ or m? -^ —rri^, or both). We anticipate that, under the approximations 
discussed below, the oscillation physics in such alternative 2+2 spectra is equivalent (under 
appropriate changes of variables) to that in the spectrum of Fig. 1. Therefore, our reference 
choice in Fig. 1 does not represent a limitation, as far as 2+2 spectra are considered.^] We 



^Although our work is focussed on 2+2 spectra, we will also briefly comments upon alternative 
3+1 spectra (a triplet plus a loner state) in Sec. Vll. In fact, the standard arguments against 3+1 



models (see, e.g., |||]) appear now to be somewhat less compelling [^lE] in light of the most recent 
LSND data M. 



also anticipate that, under the same approximations, CP violation effects turn out to be 
unobservably small in current experiments, so that the mixing matrix U can be effectively 
taken as real, 



U^ c^U^ 



(14) 



Further assumptions about the mixing matrix U can be adopted by comparing positive 
disappearance results at the mass scale Sm^ (solar z/g —>■ Ve) and m? (atmospheric v^ —>■ u^) 
with negative disappearance results at the mass scale M^ in short-baseline accelerators 



no z/„ 



i^„ 



and reactors fno Up 



[see, e.g. 



Such results imply that the z/^ 



component of the atmospheric doublet (z/3, z/4) must be large to explain the atmospheric z/^ 
anomaly, but must be small in the other doublet (z/i, 2/2), in order to avoid large M^-driven 
^^l ~^ ^fj. oscillations at short-baseline accelerators. Analogously, the u^ component of the 
solar doublet [ui, V2) must be large to explain the solar Ve deficit, but it must be small in 
the other doublet (z/3, z/4), in order to avoid large M^-driven z/g -^ z/g oscillations at reactors. 
Therefore, one can take 



U,^ + U, 






^l2 



u. 



e4 







(15) 
(16) 



Some remarks on the above two assumptions are in order. In 4z/ models embedding 
LSND, Eqs. ([TsD and (16) must be slightly violated — at the few permill level — in order to 
get a nonzero value for the small (few permill) LSND u^ -^ v^ oscillation amplitude, given 

A few permill violation makes 



LSND 



4|t/^3?7, 



e3 



+ t/„4f/, 



At4'^e4| 



4|f/^it/, 



el 



+ ^u2t/, 



At2'^e2| 



by A 

little difference in the analysis of current solar and atmospheric neutrino data, and will be 
neglected in the following. One cannot exclude, however, more general scenarios in which 
^ui + ^u2 is taken close to the weakest upper limits allowed by accelerators (of the order 
of ten percent [^]) and can thus play a nonnegligible phenomenological role 0,[16|. Such 
scenarios are beyond the scope of the present work. 

The simplifying assumptions (11^) and (|16|) lead to a specific texture for the matrix f/. 



V 



(UeX f/e2 \ 

^sl f^s2 f^s3 f^s4 

f/^3 t/^4 

Vf/rl Ur2 Ur3 Ur4 J 



(17) 



with the following implications. The state z/^ (z/g) must be a linear combination of z/3 and 
z/4 (z/i and 1^2), 



Ve = +Cujl^l + S^l'2 



(18) 
(19) 



where c = cos, s = sin, and we have introduced two mixing angles ip and uj (ranging in 
[0,7r/2]). The corresponding orthogonal combinations z/+ and z/_, defined as 



z/_ 



-s^z/3 + c^z/4 , 

-S^Z/i + C^Z/2 , 



(20) 
(21) 



must be a linear combination of the remaining flavor states v^- and Vs 









(22) 
(23) 



where ^ is the third (and last) mixing angle needed to parametrize a unitary matrix U of 
the form (p!7|). In terms of (oj,-?/',^), the matrix U reads explicitly 



U 



( c^ Su \ 

c^ s^ 



By defining the rotation matrices 



U, 



(I \ 

c^ s^ 

10 

V -s^ c^ y 



(24) 



(25) 



U, 



and 



f/. 



/I \ 
10 

c^ s^ 

V -s^ c^ y 



^ Co; S^ ^^ 

-Su c^ 

10 

V ly 



the mixing matrix U can also be written as 

where the commutation property [U^, U^] = has been used: 

Finally, it is useful to introduce a new flavor basis z/^ defined as 



jy„ 



{Ve,V^,V^,V^ 



related to the standard flavor basis by a ^-rotation, 



(26) 



(27) 



(28) 



(29) 



(30) 



^; = f/f ^o 



(31) 



In such basis, the neutrino evolution Hamiltonian reads 

n' = ^ U^^M f/J, + V2Gf UJ diag (iV„ iiV„, 0, o) U^ . (32) 

In the next subsection, it will be shown how the above Av Hamiltonian can be effectively 
decoupled in two (solar and atmospheric) 2v sub-Hamiltonians. 

Effective Hamiltonians for mixed active-sterile oscillations in matter have been recently 
discussed also in Refs. [|],|16|- We can make contact with the (pij notation of through 
the identifications oj = 0i2, ip = 034, and C, = 024, valid under the approximations 023 — 
[equivalent to Eq. (|1^)] and 0i3 ^ 0i4 ^ [equivalent to Eq. (0)]. Analogously, we can 
make contact with the 6ij notation of [0 through the identifications ^ = 6'34 and ip = 624, 
valid under the approximation 6*23 — [equivalent to Eq. (|T3])]. 

C. Approximations related to atmospheric i^'s 

As far as atmospheric neutrinos are concerned, we neglect in first approximation the 
small (z/i, z/2) mass square difference Sm'^, and take the inter-doublet mass square difference 
M^ very large. The matrix A4 in Eq. (H) for atmospheric neutrinos can then be written as 

Ma ^ diag (-M^, -M\ 0, m^) . (33) 

By inserting such A4a in Eq- (|3^), one obtains that: i) The angle u is rotated away; 
and ii) as M^ -^ cxd, the states (ue,^-) [or, equivalently, (z/i,z/2)] oscillate with very high 
frequency (oc M^), and effectively decouple from the states (2/^,//+). Therefore, only the 
evolution of (z/^, z/+) is relevant, 

. d ' ' 



I 



v^, \ „„ I V, 



nU r ' (34) 



dx \y+ ) -^ \u 

the Hamiltonian being given by 

where we have subtracted an inessential term m^/AE from the diagonal elements of 7i^. 

Equation (^) for atmospheric neutrinos is equivalent to the familiar one describing the 
i^fi -^ i^s Hamiltonian, provided that the usual mass and mixing parameters are identified 
with m? and ■0, and that the neutron density iV„ is multiplied by a factor si. For increasing 
values of s^, one gets a smooth interpolation from pure v^ -^ v^ oscillations (s| = 0) to pure 
^^ -^ ^s oscillations (s? = 1), passing through mixed active-sterile transitions (0 < s? < 1). 

Equation (pSj), which has been derived for the reference 4z/ spectrum of Fig. 1, trivially 
holds also for for spectra with negative bvr? or negative M^ (being independent of such mass 
parameters). Concerning the case of negative rn? ^ one can easily prove (by swapping the 
states 1^3 and z/4), that it is equivalent to the case of positive m,^, modulo the replacement 
s^l, — ^ c^ (and z/+ -^ — '^+), which corresponds to swap the first two octants of 0), and thus 
does not need a separate treatment (as far as -0 G [0,7r/2]). It follows that the sign of rv? 
is irrelevant at the octant boundary (-0 = vr/4), namely, at maximal atmospheric v mixing. 
For atmospheric antineutrinos, the same considerations as for neutrinos apply, provided that 
the neutron density is taken with opposite sign. 



D. Approximations related to solar z^'s 

As far as solar neutrinos are concerned, the atmospheric neutrino doublet (z/3, 1/4) can 
be seen as a couple of "far" mass eigenstates, separated from (z/i,z/2) by a large M^ gap. 
At zeroth order in 5m? /m? and in m? /M'^, the matrix M. in Eq. (§) for solar neutrinos can 
then be written as 

Ms ^ diag (0, 5m^, M^, M^) . (36) 

By inserting such Ais in Eq. (^2]), one obtains that: i) The angle ip is rotated away; 
and ii) as M^ -^ 00, the states (z/^, i/+) [or, equivalently, (1/3, i/4)] oscillate with very high 
frequency (oc M^), and effectively decouple from the states (z/e;^^-)- Therefore, only the 



evolution of (z/g, z/_) is relevant. 




n's r;] , (37) 



z/_ 



the Hamiltonian being given by 



Wi = ^f?" ^-UvIcJ'V'i'^^'"), (38) 



where we have subtracted an inessential term [dm? /AE + ^j2Gp \ ci A^„) from the diagonal 
elements of Ti'^. 

Equation (^) for solar neutrinos is equivalent to the familiar one describing two-family 
oscillations, provided that the usual mass and mixing parameters are identified with bvn? 
and cj, and that the background fermion density is taken as N ^ = A'g — |c|A^„. For increasing 
values of s\^^ one gets a smooth interpolation between pure v^. — > v^ oscillations {s\^ = and 
Nf = Ne — |iV„) and pure z/g — ^ i^t oscillations (s| = 1 and Nf = N^), passing through 
mixed active-sterile transitions (0 < s| < 1). 

Equation (pSD, which has been derived for the 4z/ reference spectrum of Fig. 1, trivially 
holds also for spectra with negative m^ or negative M^ (being independent of such mass 
parameters). Concerning the case of negative 6ir?, one can easily prove (by swapping the 
states z/i and z/2), that it is equivalent to the case of positive 6m?, modulo the replacement 
Suj —>■ c^ (and z/_ —>■ —i^-), which correspond to swap the first two octants of u, and thus 
does not need a separate treatment (as far as cj G [0, vr/2]). It follows that the sign of 6m? 
is irrelevant at the octant boundary [uj = vr/4), namely, at maximal solar z/ mixing. 



E. Summary and remarks about the 4z^ framework 

We have introduced a 2+2 four-neutrino scenario, based on a few phenomenological 
approximations, which embeds at the same time both active and sterile oscillations of at- 
mospheric and solar neutrinos. The final states z/+ and z/_ for atmospheric and solar flavor 
transitions (z/^ -^ z/+ and u^ -^ z/_, respectively) represent orthogonal combinations of Ur 
and z/g through a mixing angle ^ [Eqs. (|2^ ) and (pS])]. 



The angle ^ governs the relative amount of Us and Ur in the states z/±, and thus also 
the amplitude of matter effects. When ^ ^ 0, the state u^ (oscillation partner of u^) is 
dominantly a z/^, while the state z/_ (oscillation partner of z/g) is dominantly a Ug] and vice 
versa for ^ -^ n/2. For ^ = 7r/4, both solar and atmospheric neutrino oscillations are 
simultaneously and democratically distributed into the active and sterile channels. 

Atmospheric neutrino oscillations are governed by an effective 2z/ hamiltonian with vac- 
uum mass-mixing parameters (m^, ip), and with an effective fermion density in matter given 
by Nf = ^s'^Nn [Eq. (pS])]. The case of negative m^ is equivalent to swap the first two 
octants of V', and thus it does not need a separate treatment. 

Solar neutrino oscillations are governed by an effective 2i/ hamiltonian with vacuum 
mass-mixing parameters {6m'^,uj), and with an effective fermion density in matter given by 
Nf = Ne — ^c^Nn [Eq. (|38|)]. The case of negative 5m^ is equivalent to swap the first two 
octants of uj, and thus it does not need a separate treatment. 

The sign of M^ (namely, the occurrence of a solar u doublet heavier or lighter than the 
atmospheric u doublet) is irrelevant under our approximations, as far as current oscillation 
experiments are concerned. Such sign can instead be relevant in other contexts, such as in 
supernovae (see, e.g., |21|) in big-bang nucleosynthesis (BBN) (see, e.g., p2| , |23[| and refs. 
therein) and, for Majorana neutrinos, in neutrinoless double beta decay (Oz/2/?) (see, e.g., 
|p^^6| and refs. therein). In particular, BBN data can probe new (sterile) neutrino species 
in addition to the three active neutrinos, and thus constrain Au scenarios. However, at 



present, such constraints do not appear to be stable enough |2^ to really exclude specific Au 
models with great confidence, especially if one elaborates upon a possible tension between 
the most recent BBN and cosmic microwave background data ||28|. Moreover, subtle effects 



of sterile neutrino oscillations in the early universe [2^, which are still subject to intensive 



studies (see [Q and refs. therein) might provide ways to evade or mitigate standard BBN 
bounds. Constraints from Oi/2/5 decay can also be evaded by assuming Dirac neutrinos. For 
such reasons, we postpone the discussion of (potentially important but presently uncertain) 
additional BBN and Oz/2/5 bounds to a future work, while in this paper we prefer to focus 
on a homogeneous set of neutrino oscillation data (SK atmospheric). 

Finally, we point out that the reduction from the general 4z/ oscillation Hamiltonian 
(^) to effective 2z/ forms for both atmospheric and solar neutrinos [Eqs. (|35| ) and (^8|), 
respectively] justifies a posteriori our neglect of CP violation phases (unobservable in 2z/ 
scenarios), and thus the initial position (p!^ . An alternative proof of the vanishing of 
(observable) CP violation effects in our framework can also be obtained by applying our 
approximations ([T3|), (0) and (|^) to the 4z/ Jarlskog invariants explictly worked out in 
PI. 



III. GRAPHICAL REPRESENTATIONS 

In our Au framework, it turns out that the mixing parameter spaces for atmospheric and 
solar neutrinos can be usefully represented in triangular plots, embedding unitarity relations 
of the kind UI + U2+UI = 1 through the heights projected within a triangle with equal sides 



(and unit total height). Such unitarity relations hold for the four columns of the mixing 
matrix U characterizing our framework [Eq. (|l^)], and can be separated into two constraints 
on the atmospheric (1^3,1^4) doublet, 

U'^, + Ul, + Ul = 1 , (39) 

Ul, + Ul, + Ul = 1 , (40) 

and two constraints on the solar (z^i, U2) doublet, 

Uj, + U\ + Ul = 1 , (41) 

UI2 + Ul^ + Ul = 1 . (42) 

As we will see below, it is sufficient to implement just one unitarity relation in one triangle 
plot for each doublet, the other relation being a consequence. We choose to implement (^OD 
for atmospheric neutrinos and (|4^) for solar neutrinos. 

Figure 2 represents the "atmospheric neutrino triangle," embedding the unitarity con- 
straint (^0[) , which is related to the flavor composition of z/4, 

Z/4 = U^i Uf, + UrA ^T + UsA Us (43) 

= S^Uf, + C^C^ Vr + c^s^ Us (44) 

= s^Uf, + c^ u+ . (45) 

In the triangle plot, the upper, lower left and lower right corner are identified with u^, z/,-, 
and Ug, respectively. The heights projected from a generic point z/4 inside the triangle onto 
the lower, right, and left side represent the elements U'^^, U^^, and 11"^^, respectively. When 
ua coincides with one of the corners (mass eigenstate = fiavor eigenstate), no oscillation oc- 
curs. Generic inner points describe mixed (active+sterile) atmospheric neutrino oscillations, 
smoothly interpolating from pure z/^ -^ v^ (left side) to pure z/^ -^ Ug (right side). In the 
upper triangle plot of Fig. 2, the (s^, s'^) parametrization is also charted. The lower triangle 
plot shows explicitly that z/4 is a linear combination of z/^ and z/+, with z/+ confined to the 
lower side (being a linear combination of z/,- and z/^). Notice that, once the point z/4 = (s?, si) 
is fixed, the position of z/3 is also determined in the same plot at the coordinates {sj,c^), 
and there is no need for a separate triangle embedding the constraint (|39|) . The points z/3 
and z/4 are then symmetrically placed onto the line joining the states z/^ and z/4. (of which 
they are orthogonal combinations). 

Figure 3 represents the "solar neutrino triangle," embedding the unitarity constraint 
(1^), which is related to the fiavor composition of z/2, 

Z/2 = Ue2 Ve + ^t2 ^t + V ^2 ^s (46) 

= S^ Z/e - C^Sg Z/^ + C^Sg Vs (47) 

= So; Z/e + Co; Z/_ . (48) 



^Triangle plots have been already introduced and discussed in detail in the context of 3z/ mixing, 
see Hg. 
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In the triangle plot, the upper, lower left and lower right corner are identified with i/g, Us, 
and z/t-, respectively. The heights projected from a generic point 1/2 inside the triangle onto 
the lower, right, and left side represent the elements 17^2^ U'^ai ^^^ L^^4, respectively. When U2 
coincides with one of the corners (mass eigenstate = fiavor eigenstate), no oscillation occur. 
Generic inner points describe mixed (active+sterile) solar neutrino oscillations, smoothly 
interpolating from pure v^ —^ s (left side) to pure v^ -^ v^ (right side). In the upper triangle 
plot of Fig 3, the (s^, s^) parametrization is also charted. The lower triangle plot shows 
explicitly that z/2 is a linear combination of z/g and v_^ with v_ confined to the lower side 
(being a linear combination of z/^ and v^). Notice that, once the point v^ = {sj, s^) is fixed, 
the position of z/i is also determined in the same plot with coordinates (s|, c^), and there is 
no need for a separate triangle embedding the constraint (^l]). The points z/i and 1^2 are then 



symmetrically placed onto the line joining the states z/g and z/_ (of which they are orthogonal 
combinations) . 

Finally, Fig. 4 shows the link between the solar and atmospheric parameter spaces, 
which follows from the fact that the variable sin^ ^ must represent the same abscissa in both 
triangles of Figs. 2 and 3. By putting the two triangles on top of each other, it follows that 
the Z/+ and z/_ points must be placed on the same vertical line, and that couple of points 
inside the atmospheric and solar triangles must be placed on iso-s| lines (slanted dashed 
lines in Fig. 4). It will be shown that such "common ^" constraint will significantly reduce 
the solutions that would be separately allowed by solar data only or by atmospheric data 
only. 

IV. SUPER-KAMIOKANDE ATMOSPHERIC u DATA AND EXPECTATIONS 

In the previous sections we have set the theoretical framework, as well as the graphical 
representations, for our analysis of 4z/ oscillation solutions to the atmospheric u anomaly. 
In this section we present the data used in the fit, and make a first discussion (to be refined 
later) of the relative importance of different data samples in constraining the mixing of 
atmospheric z/^ with Ug. 

Figure 5 shows the Super-Kamiokande data used in our x^ analysis, as graphically re- 
duced from 0, corresponding to a detector exposure of 70.5 kTy. The reported errors are 
statistical only (±lcr). Systematic errors are treated as in ^^. The data are shown as binned 



distributions of e-like and /i-like event rates R, normalized to no-oscillation expectations Rq 
in each bin, in terms of the lepton zenith angle i). The distributions refer to sub-GeV (SG) 
leptons, multi-GeV (MG) leptons, upward stopping (US) muons, and upward through-going 
(UT) muons, according to the SK terminology |)7||52|. The total number of data points used 



in the fit is 55. The SK data on neutral-current enriched samples |q], which are also useful 



in testing the presence of a Ug component |3^, are not included in the present analysis. 



Although our main results are based on accurate calculations [^ of the full zenith 
distributions in Fig. 5, it is useful to gain some prior understanding of the physics by 
means of (partly) averaged quantities such as the popular fi/e double fiavor ratio for SG 
and MG events, the up-down asymmetry of MG events, and the vertical-to-horizontal ratio 
of UT events. We follow the SK conventions for such quantities D^J^, namely, data with 
I cos-i^l < 0.2 are removed in the up-down asymmetry, and the separation between "vertical" 
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and "horizontal" through-going muons is conventionally taken at cos^ = —0.4. 

Figure 6 shows the double ratio R^/e of MG events vs SG events. The corresponding 
SK data are shown as a cross (with statistical and total la error bars). The theoretical 
expectations at maximal mixing are shown as a solid line for pure i/^ — » Vr oscillations 
(s| = in our notation) and as a dashed line for pure z/^ —* v^ oscillations (s| = 1). Such 
lines connect the points calculated at ±m^ = n x 10*" eV^, with m = —5, —4, . . . , —1, 
and n = 1, 2, . . . , 9 (the sign of m^ is irrelevant at maximal mixing, s| = 1/2). Cases with 
< s| < 1 (not shown) would give lines intermediate between the solid and the dashed ones. 
The "trajectories" of the theoretical points start at R^/e — 1 for small m^ (no-oscillation 
limit), and end at -R^/e ~ 1/2 for large m^ (averaged oscillations). 

From Fig. 6 we learn that the double ratio is not a good variable to discriminate between 
J^^l — *■ ^T and z/^ -^ Ug oscillations, since the separation of the two corresponding curves is 
much smaller than the experimental error. In fact, the matter effects which should distin- 
guish the Vg channel are suppressed both by the relatively low energy of SG and MG events 
and by integration over the zenith angle. On the other hand, the combined SG and MG 
information on R^/e in Fig. 6 is sensitive to m^, and favors the range ~ [10^^, 10^^] eV^, 
with ~ 4 X 10~^ eV^ being the value closest to the data for both oscillation channels. We 
conclude that the SG and MG yu/e double ratios give valuable information on m? , rather 
than on the flavor of the z/^ oscillation partner. 

A better sensitivity to the sterile component of the t'^ oscillation partner can be gained by 
improving the zenith angle information and by probing higher energies, where matter effects 



are larger |^. Figure 7 shows the ratio of vertical to horizontal (V/H) upward through- 
going muons vs the up-down asymmetry of multi-GeV muons {U — D /U + D). As in Fig. 6, 
solid (dashed) lines refer to u^ -^ v^ {u^ -^ Ug) oscillations at maximal mixing, s^ = 1/2. 
In addition. Fig. 7 shows, as a dotted line, the case of mixed active-sterile oscillations with 
equal amplitude [z/+ = (z/^ + h'r)/V^, namely, s| = 1/2]. The trajectories of the theoretical 
points start and end at small values of the up-down asymmetry, as expected in the limits 
m^ — > (no oscillation) and m^ —>■ cxo (averaged oscillations), after passing through negative 
values (corresponding to a suppression of the upgoing muon rate). The V/H ratio should 
also take back its no-oscillation value as m^ -^ oo, although the highest values reported 
in Fig. 7 {m'^ = 10~^ eV^) is not "asymptotic" enough to show this behavior. The V/H 
variations with m^ are rather small in the pure u^ -^ z/^ case, when compared with the 
corresponding variations in the pure z/^ — >• u^. case, as a consequence of the strong damping 
of the sterile transition amplitude in matter at high energies [^ . Intermediate variations of 
V/H occur for the mixed active-(-sterile case (dotted line); therefore, the V/H asymmetry 
is a sensitive probe of s| through the effects of the matter term ^Gp s| A^n/2. 

The recent SK data 0, shown in Fig. 7 as a cross with error bars (statistical and total 
at Icr) clearly favor the pure z/^ — > v^ case (solid line), as compared with the pure z/^ — *■ Vg 
case (dashed line). By themselves, such data cannot exclude an intermediate situation with 
large mixing of v^ with v^ and Vg (dotted line). In such case, however, the favored values of 
m? are rather large [~ O(10~^] eV^), in contrast with the information coming from the /i/e 
ratio in Fig. 6. Therefore, the interplay of different pieces of data seems to indicate clearly 
an upper bound on the Vg component of z/+ (namely, on s|). 
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V. RESULTS OF THE ATMOSPHERIC u ANALYSIS 

In this section we present the results of our fit to the atmospheric neutrino data of Fig. 5 
(55 data points), as obtained by calculating the SK zenith distributions for unconstrained 
values of the three relevant parameters (m^, s|, s^). We discuss first the bounds on the mass 
parameter m^, and then those on the mixing parameters (s|, sV). 

Figure 8 shows the value of the global x^ as a function of m^ (in linear scale), for 
four z/^ — i> Z/+ oscillation cases, corresponding to: {i) unconstrained z/+ (thick solid line); 
(u) Z/+ = Vr (thin solid line); {Hi) z/+ = Vs (dashed line); and {iv) z/+ = (z/^ + z/t-)/\/2 
(dotted line). In terms of the mixing angle ^, such cases correspond, respectively, to: {i) 
unconstrained s|; {ii) s| = 0; (m) s| = 1; and {iv) s| = 1/2. The parameter s^ is left 
free in all four cases. The global minimum for the unconstrained case (Xmin ~ 47.4) is 
reached at (m^/eV^, s|, s^) = (3.2 x 10"'^, 0.18, 0.51). The slight preference for a nonzero 
value of sj at the point of minimum x^ is intriguing, but not statistically significant, since 
the fit with s| = (pure v^ —>■ Vj. oscillations, thin solid line) is only slightly worse than the 
unconstrained fit (thick solid line). 

Concerning the statistical interpretation of the results shown in Fig. 8, some remarks are 
in order. In general, one can use the x^ function for two complementary purposes [^]: (1) 
to estimate the best-fit values of {m?, sj, s\) and their uncertainties {parameter estimation)] 
or (2) to test the goodness-of-fit with some prior assumptions about (m^, s|, s^) {hypothesis 
test). In the first approach, the values of (m^,s|,s3,) minimizing the x^ function are taken 
as best estimates of the (unknown) true values of such parameters. The corresponding 
uncertainties are then obtained by allowing a shift in the x^ function (Ax^ = X^ ~ Xmin); 
with corresponding confidence levels determined by the number of free parameters in the 
model {Ndf = 3 in our case). From this point of view, the case of pure z/^ — > Us oscillations, 
corresponding to Ax'^/N^p = 14.7/3 in Fig. 8, is reached only at 99.8% C.L., and is thus 
strongly disfavored. In the second approach, one checks how well a specific hypothesis 
about the parameters fits the 55 SK data points, independently of the existence of a global 
X^ minimum at some point in the parameter space. For instance, if order to test the a priori 
hypotheses of pure z/^ -^ Ur or pure z/^ — > Ug oscillations, one should compare the absolute 
X^ value with A^df = 55 — 2, the two free variables being m^ and s^ (s| being fixed at or 1, 
respectively). From Fig. 8, we then get that the probability of having a worse x^ is 63% for 
s| = (pure u^ -^ z/^ oscillations, x^/Ndf = 49/53) and 18% for s| = 1 (xV^Vdf = 49/53 
(pure Ufj, -^ Vs oscillations, x^/^df = 62/53). Therefore, although the pure sterile case is 
clearly disfavored in comparison with the pure active case, it is not ruled out a priori in 
our oscillation analysis. Having clarified such statistical issues, we proceed to estimate the 
mass-mixing parameters {rn^ , s'j., sV) by using the first approach, based on Ax^- 

In Fig. 8, all four cases show a preference for rn^ ~ 3 x 10~^ eV^, which thus represents 
a very stable indication coming from the SK data, independently of the amplitude of the Us 
component. For the most general case (unconstrained v+) we obtain 



m 



2 ^ Q 0+2.3 ,, 1 n-3 „-\/2 



3.2+f:^ X 10-^ eV' (90% C.L.) , (49) 

by taking Ax^ = 6.25 for iVoF = 3. 

Figure 8 also gives gives indications about the parameter s"^^. The case s| = (thin 
solid line) gives results comparable to the best fit for unconstrained s"^^ (thick solid line). 

13 



The case of a fifty-fifty admixture of Ug and Ut in z/+ {sj = 1/2, dotted line) provides a fit 
that, although worse than for pure z/^ — ;> i/,-, is still acceptable (its minimum x^ being only 
four units above the absolute minimum). However, the quality of the fit worsen rapidly 
as Z/+ gets closer to Ug, the worst case being reached for pure u^ -^ Ug oscillations (dashed 
line). Therefore, we expect to get an upper bound on sj from the analysis of the mixing 
parameters, which can be clearly represented in the triangle plot. 

Figures 9 and 10 show the results of our analysis in the atmospheric triangle plot, for 
10 representative (decreasing) values of m^. The first three columns of triangles refer to the 
separate fits to sub-GeV electrons and muons (10+10 bins), multi-GeV electrons and muons 
(10+10 bins), and upward stopping and through-going muons (5+10 bins), while the fourth 
column refers to the total SK data sample (55 bins). For each column, we find the minimum 
X^ from the fit to the corresponding data sample, and then present sections of the allowed 
volume at fixed values of rn^ for Ax^ = 6.25 (90% C.L., sohd lines) and Ax^ = 11.36 (99% 
C.L., dotted lines). It can be seen that the low-energy SG data are basically insensitive to 
s|, since they are consistent both with z/^ —^ v^ oscillations (left side) and with z/^ -^ Vg 
oscillations (right side), as well as with any intermediate combination of the two oscillation 
channels. High-energy upgoing /i data are instead much more sensitive to the Vg component 
through matter effects, which increase both with energy and with s|, and tend to suppress 



the /i deficit [^. Large matter effects appear to be disfavored by the upgoing muon data, 
the rightmost part of the triangle being excluded at 90% C.L. The multi-GeV data cover 
an intermediate energy range and are not as constraining as the upgoing muons; however, 
they also show a tendency to disfavor a large Vg component. This tendency is strengthened 
in the global combination of data (fourth column), leading to the upper bound 

s| < 0.67 (90% C.L.) , (50) 

which clearly disfavors pure or quasi-pure v^ -^ Ug oscillations. In particular, from our 
Figs. 9 and 10 we derive that pure u^ — > Ug oscillations (s| = 1) are excluded at > 99% C.L., 
consistently with the SK analysis 0,|]. 

However, the bound (|50|) still allows intermediate cases of combined active+sterile os- 
cillations with a sizable Ug component. For instance, one cannot exclude that the sterile 
neutrino channel may have the same amplitude as the active one (s| = 0.5), or may even be 
dominant (e.g., s? = 0.6, corresponding to z/_|_ ~ 0.77 z/^ + 0.65 z/,-). Furthermore, for m^ in 
its upper range (e.g., m^ = 5 x 10~^ eV^ in Fig. 9), a subdominant Ug component actually 
helps in the global fit to SK data, leading to a 90% C.L. allowed region which does not touch 
the left side of the triangle (corresponding to the pure u^ -^ z/,- channel). 

Figures 9 and 10 also provide bounds on the other mixing parameter s'i, which governs the 
relative amount of z/^ and u^ on the atmospheric neutrino states z/3 and z/4 [see Eqs. 



and Fig. 2], and thus the overall amplitude of z/^ -^ z/+ oscillations. From Figs. 9 and 10 we 
derive 

4^0.51 + 0.17 (90% C.L.) , (51) 

which favors z/^ — >■ z/+ oscillations with nearly maximal amplitude, as one expects from the 
observation of nearly maximal average suppression (~ 50%) of upgoing MG muons. There 
is a slight asymmetry of the allowed regions with respect to s^ = 1/2, which reflects the 
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relative change of sign of vacuum and matter terms in Eq. (^) when passing from the first 
to the second octant. 

Figures 8, 9, and 10 demonstrate that, on the basis of present SK data on the zenith 
distributions of leptons induced by atmospheric neutrinos, pure z/^ — > Ug oscillations are 
strongly disfavored, but one cannot exclude mixed active-sterile z/^ -^ u^ oscillations with 
I/+ = c^Ur + s^Ug, provided that the partial amplitude of the sterile channel is ^ 67% 
[Eq. (ISDD]. Since such bound is dominated by high-energy muons, it is important to under- 
stand how such data can improve the sensitivity to s| in the future. 

Figure 11 shows the zenith distributions of upward through-going muons at maximal 
(z/^, z/_,_) mixing (s^ = 1/2) for z/+ = z/^ (left panel), z/+ = (z/,- + Us)/V2 (middle panel), 
and Z/+ = Us (right panel). The solid (dashed) histograms refer to ±?ti^ = 3 x 10~^ eV^ 
(5 X 10^^ eV^), the sign of m"^ being irrelevant at maximal mixing. From left to right, sj 
increases and matter effects also grow with it, producing an increasing modulation of the 
zenith distributions. The modulation, which is very strong for the oscillation probability 
at fixed energies (not shown), is largely smeared out in Fig. 11 by integration over the 
energy spectrum and, to a lesser extent, by the presence of both z/ and z7 (with opposite 
matter terms) in the atmospheric neutrino flux. Therefore, the modulation pattern of the 
UT/i distribution might be more clearly observed in future atmospheric z/ experiments ||^ 



aiming at improved energy reconstruction and /i+//i~ separation. Concerning the present 
SK data (Fig. 11), the large scatter of central values and the still large error bars prevent 
a clear discrimination of the theoretical distributions in the left and middle panels, which 
anyway appear to provide a better fit than in right panel. A reduction by a factor ^ 2 in the 
uncertainties of the SK UT/i sample appears necessary for a decisive discrimination among 
the three cases shown in Fig. 11. In conclusion, we can still learn a lot on the Ug component 
of Z/+ as we get more high-energy muon data in SK, as well as in MACRO, and possibly 
from future atmospheric p experiments. 

VI. LINKING ATMOSPHERIC AND SOLAR NEUTRINOS 

Solar neutrino oscillations with both active and sterile states have been recently investi- 



gated in ||T3|Jl^ , whose results can be interpreted, in our framework, through the identifica- 



tion CI3C24 = 1 — s|, 6*12 = uj, and Am'^i = 5m?. In fT^JT^ it has been shown that there is a 
continuous set of solutions to the solar neutrino problem for unconstrained s|, ranging from 
pure "active oscillation solutions" (s| = 1) to pure "sterile oscillation solutions" (s| = 0). 
In particular, for si = 1 all the usual solutions to the solar neutrino problem are present, 
either at small uj mixing (matter-enhanced solution) or at large u mixing (matter-enhanced 
solution at low or high 5m^, or vacuum solutions). For decreasing values of s^^, all the so- 
lutions tend to shrink; however, the large mixing (LM) ones (vacuum or matter-enhanced) 
eventually disappear for s| < 0.3-0.4, while the small mixing (SM) solution still survives 
(at ~ 99% C.L.) even for si = 0. Such behavior is qualitatively reported in Fig. 12, through 
the triangular representation of the parameter space. The plot is only qualitative, since the 
detailed behavior of any solution depends both on the confidence level adopted and on the 
specific value chosen for 6171'^, whose complete scan is beyond the scope of this paper. 

In Fig. 12, the allowed region at large mixing (s^ close to 1/2), which has the maximum 
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likelihood at s| = 1 (pure active oscillations, right side of the triangle), gradually shrinks 
when sj decreases and eventually disappears. The allowed region at small mixing (s^ close 
to 0) covers instead the whole s| range. 

As discussed in Sees. II and III, the solar and atmospheric solutions are coupled through 
the common angle ^. Such coupling can be visualized by putting the solar and atmospheric 
triangles on top of each other, as shown in Fig. 4. The results are qualitatively different for 
the SM and LM solutions to the solar neutrino problem, as shown in Figs. 13 and 14. 

Figure 13 shows, in grey, the regions globally allowed by the combination of solar and 
atmospheric z/ data, assuming the SM solution to the solar u problem. The fact that atmo- 
spheric u data put an upper bound on s| (as quantitatively discussed through Figs. 9 and 
10) implies that also the fraction of the SM solar z/ solution at large s| is excluded. However, 
there is still a continuous range of global solutions for s| 7^ 0, having a sizable Us component 
for both solar and atmospheric z/ oscillations. Therefore atmospheric u data "cut" a frac- 
tion of the SM solution to the solar z/ problem, but the combination still allows solutions 
involving active+sterile neutrino oscillations for both atmospheric and solar neutrinos. 

The opposite happens for the LM solution to the solar neutrino problem. In this case, 
solar neutrino data cut the "small s|" part of the atmospheric allowed region, as shown 
in Fig. 14, and the globally allowed solutions (grey inner regions) represent genuine ac- 
tive-|-sterile oscillations, which never reduce to the pure active or pure sterile subcases (the 
sides of the triangles). In particular, the points denoted by stars in Fig. 14 represent a 
peculiar solution (which might be dubbed "fourfold maximal mixing"), corresponding to 
maximal mixing of z/± with z/^.t (s| = 1/2), maximal mixing of z/3,4 with z/^^+ (s^ = 1/2), 
and maximal mixing of z/i_2 with z/e^_ (s^ = 1/2). In this case, both solar and atmospheric 
neutrino oscillations are equally distributed among the active and sterile channels. 

In conclusion, it turns out that, at least qualitatively, 4z/ scenarios with mixed ac- 
tive+sterile transitions for both solar and atmospheric neutrinos are allowed by the com- 
bination of world oscillation data, with different mixing patterns dictated by the solution 
assumed for the solar neutrino problem (either small or large mixing). Consistency between 
atmospheric and solar neutrino data in 4z/ models has also been recently found in an inde- 
pendent and somewhat different analysis []16[. Therefore, the hypothesis of a fourth, sterile 



neutrino is still compatible with present data and needs further theoretical and experimental 
investigations. Such results clearly demonstrate that the two-flavor approximations often 
used to analyze solar or atmospheric neutrinos data (in terms of either pure active or pure 
sterile oscillations) are, in general, not appropriate to explore the full parameter space al- 
lowed by 4z/ models. Our 4z/ formalism and triangular representations can also be usefully 
applied to explore the z/^ -^ Ug discovery potential in long baseline accelerator experiments, 
as shown for the MINOS experiment in 



VII. REMARKS ON 3-M MASS SPECTRA 

The standard arguments disfavoring 3+1 spectra (triplet+singlet) with respect to 2+2 
spectra (two doublets) in the interpretation of the current 4z/ oscillation phenomenology |l| 
appear to be somewhat weakened ||T^,[TP| in light of the most recent LSND data 0, and 



should be perhaps revisited quantitatively. Here we make just a few qualitative comments 
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on the mixing pattern that is hkely to emerge in 3+1 models characterized by a triplet 
(z/i, z/2, z/3) (mainly responsible for solar and atmospheric oscillations) and by a "loner" state 
z/4 (mainly responsible for LSND oscillations). 

The absence (or weak occurrence, as in LSND) of oscillations in all channels probed 
by short-baseline reactor and accelerator experiments is consistent with 1/4 being close to a 
flavor eigenstate. The dominant flavor of z/4 cannot be Ve i^^/j,), otherwise the disappearance 
pattern of solar z/g (atmospheric u^) would remain unexplained. Moreover, if it were z/4 ~ z/,-, 
then the triplet (z/i, z/2, z/3) would approximately be a linear combination of (z/^, z/g, z^s), and 
thus the atmospheric z/^ oscillation partner should be either p^ (excluded) or u^ (disfavored 
at ^ 99% CL.) or a linear combination of z/g and z/^ (presumably disfavored). Therefore, the 
favored option for 3+1 spectra seems to imply z/4 ~ z/^, with (at least) a small subdominant 
component of z/^ and z/g in z/4, necessary to drive LSND oscillations (see also [p!^,^ ) f| 

Assuming z/4 ~ z/^ in 3+1 models, the triplet states (z/i, z/2, z/3) would then approximately 
be linear combinations of (z/g, z/^, z^r), and the triplet phenomenology of solar and atmospheric 
neutrinos would practically become equivalent to that of a standard framework with three 
active neutrinos. Such 3z/ framework has been worked out in detail in previous papers 
|3^ , |33| , ^ , where the reader can find the corresponding constraints on the mass-mixing 



parameters. In conclusion, it seems that the emerging mixing pattern for the 3+1 spectrum 
implies, up to small "LSND perturbations", an effective three-neutrino phenomenology for 
solar and atmospheric oscillations. In this case, it becomes important to study how well one 
can experimentally discriminate such 4z/ option from pure 3z/ oscillations. 

VIII. SUMMARY AND CONCLUSIONS 

In the context of 2+2 neutrino models, we have thoroughly analyzed the SK atmospheric 
data on the zenith distributions of lepton events, assuming the coexistence of both z/^ — >■ z/t- 
and z/^ — i> Us oscillations, with a smooth interpolation between such two subcases. We have 
shown that, although the data disfavor oscillations in the pure u^ -^ z/^ channel, one cannot 
exclude their presence in addition to z/^ — > z/^- oscillations. High energy muon data appear 
to be crucial to assess the relative amplitude of the active and sterile oscillation channels for 
atmospheric neutrinos. We have also shown, in a qualitative way, that mixed active+sterile 
oscillations of atmospheric neutrinos are compatible with analogous oscillations of solar 
neutrinos, with different mixing pattern emerging from different choices for the solar z/ 
solutions. Such results indicate that the solar and atmospheric neutrino phenomenology 
cannot be simply embedded in the usual two-family scenarios involving oscillations into 
either pure active or pure sterile states, and call for an extended 4z/ framework involving 
combined solutions. Further directions to improve and constrain our present results might 
involve, on the one hand, a quantitative, global analysis of world neutrino oscillation data 
and, on the other hand, the inclusion of non-oscillatory data which have not been considered 
in this work, such as those coming from astrophysics and cosmology, and from direct neutrino 
mass measurements. 



^Alternative 3+1 models with small (2/^,1/4) mixing are currently being investigated, see 
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Reference 4i/ mass spectrum 
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Fig. 1. The reference 4zv mass spectrum adopted in this work. 
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Atmospheric v parameter space 
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Fig. 2. The atmospheric v parameter space, in the triangular representation. See the text for 
details. 
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Fig. 3. The solar u parameter space, in the triangular representation. See the text for details. 
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Solor - Atmospheric link 
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Fig. 4. Graphical representation of the hnk existing between the solar and atmospheric param- 
eter spaces. 
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Super-Kamiokande (70.5 kTy) 
e, fx zenith distributions 
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Fig. 5. Super-Kaniiokande data on the zenith distributions of lepton events induced by atmo- 
spheric neutrinos. 
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yL/-/e double ratio of sub-GeV and multl-GeV events 



0.9 



0.8 



R 



MG 



0.7 



0.6 



0.5 



0.4 



7 '-^ 1 9 

^ ±m' = 10"' eV' 



1 0"' eV' 




e — B-#-B 



Theory {2u max. mixing) 

• u,, ^> jy^ 



i^M^^^s 



0.4 



0.5 



0.6 



0.7 0. 

SG 



0.9 



R 



/u,e 



Fig. 6. Double ratio of rtiuon-to-electron events for multi-GeV events (y-axis) and sub-GeV 
events (x-axis). The theoretical expectation at fixed (maximal) mixing and for running m? is shown 
as a solid curve for pure active oscillations, and as a dashed curve for pure sterile oscillations. The 
cross with error bars (statistical and total) represents the SK data. 



26 



Muon distribution asymmetries for UT and MG events 
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Fig. 7. Vertical-to-horizontal ratio of upgoing muon rates vs the up-down asymmetry of multi- 
GeV muon events. The theoretical expectation at fixed (maximal) mixing and for running m? is 
shown as a solid curve for pure active oscillations, as a dashed curve for pure sterile oscillations, 
and as a dotted curve for mixed active+sterile oscillations of equal amplitude. The cross with error 
bars represents the SK data. 
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Atmospheric v^ — > i/+ oscillation fit 
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Fig. 8. Result of the global x^ fit to the 55 data points of Fig. 5, as a function of the mass 
parameter ttt,^, assuming v^ -^ u^ oscillations for atmospheric neutrinos, and different options for 
Z/+. The mixing angle s3, is left unconstrained in all cases. The case z^+ = Vs provides the worst 
fit, but one cannot exclude a sizable Ug component of z/+. 
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m^ (eV^) Sub-GeV Multi-GeV Up-going SK Combined 
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Fig. 9. Results of the global x^ fit to atmospheric neutrino data in the triangular representation, 
for decreasing values of m^. See the text for details. 
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m^ (eV^) Sub-GeV Multi-GeV Up-going SK Combined 
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Fig. 10. As in Fig. 9, but for lower values of m^. 
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SK UJ/u distribution (1138doys) 
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Fig. 11. The zenith distribution of upward through- going muon events for three different 
oscillation scenarios, involving pure active oscillations (left panel), pure sterile oscillations (right 
panel) and mixed active+sterile oscillations (middle panel) at maximal mixing. 
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Solar V solutions (qualitative 
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Fig. 12. Qualitative pattern of \v oscillation solutions to the solar neutrino problem at small 
or large mixing, in the triangular representation of the parameter space. 
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Atm. + Solar SM (qualitative) 
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Fig. 13. Qualitative combination of the atmospheric v solution with a typical small-angle solar 
V solution. 
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Atm. + Solar LM (qualitative) 
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Fig. 14. Qualitative combination of the atmospheric i' solution with a typical large-angle solar 
i^ solution. 
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